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1 Thomas- Fermi Theory 

The Thomas-Fermi method |l|] was designed for the calculation of the electron density in 
a heavy atom, by treating the electrons as locally free. Lieb and Simon P] showed that the 
treatment is exact in the limit when the atomic number goes to infinity. Application to a 
confined Bose condensate was pioneered by Goldman, Silvera, and Legget [|], and by Oliva 

and recently reconsidered by Chou, Yang, and Yu 0. I shall describe some work on this 
subject, done in collaboration with E. Timmermans and P. Tommasini [0. 

First, let us review the original method of Thomas and Fermi. Suppose V{r) = — e$(r) 
denotes the effective potential energy of an electron in an atom, at a distance r from the 
nucleus. (See Fig.l). 
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Figure 1: Potential energy of an electron in atom. 



The condition that the electron is in a bound orbit is that 



2m 



+ Vir)<0 {!] 



^Contribution to Workshop on Bose-Einstein Condensation, Institute for Theoretical Atomic and Molec- 
ular Physics, Harvard-Smithsonian Center for Astrophysics, Cambridge, MA, (August 19-30, 1996). 
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where k is the local wave vector. Assume that all available states are occupied. Then the 
local Fermi momentum hkp{r) is given by 



2m 

which gives 

hkp = ^2me$(r) 
We know that kpi^r) is related to tke local density n(r) by 

2(47r/3)4 = n{r) 



(2) 
(3) 
(4) 



where the factor 2 comes from spin. This relates the density to the potential <l>(r). We now 
use the Poisson equation 

V2$(r) = -An[Ze5{r) - en(r)] (5) 
For r 7^ this equation is of the form 



V2$(r) +C$^/2(r) = 



(6) 



where C is a constant. One can solve this equation, and then obtain n(r). A comparison 
between Thomas-Fermi and Hartree results for Rb is sketched in Fig. 2. 
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Figure 2: Electron density in Rb: Comparison between Thomas- Fermi and Hartree approx- 
imations. 



The essence of the method is that one assumes there is a local chemical potential fJ-esi^), 
related to the true chemical potential /i by 

/Xeff(r) = /i-F(r) (7) 

In the earlier discussion, the chemical potential (Fermi energy) was taken to be zero. 
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2 Ideal Bose Gas in External Potential 

Can we apply this idea to a Bose gas? Let us first concsider an ideal Bose gas in an external 
potential, at a temperature T above the transition point. Take the potential to be harmonic: 

V{r) = ^ma;V^ (8) 
For the Thomas-Fermi idea to be applicable, we require 

Above the transition temperature, the density is related to the fugacity z — exp{—ii/kT) 
through 

n = ^gsMz) (10) 
where A = /mkT is the thermal wavelength, and 

oo 

^3/2(^)=E^ (11) 
£=1 ^ 

This suggests that in the presence of an external potential we take the local density to be 

n{r) = (^e-^^^"^)) (12) 

where /3 = 1 / kT. Integrating both sides over all space, we obtain an expression for the total 
number of particles: 

N^yJ Srg,„ (.e-^^W) (13) 

We know that g3/2{z) is bounded for < 2; < 1. So the right side is bounded. This 
forces Bose-Einstein condensation when N exceeds the bound. The number of atoms in the 
condensate A^o is given through 

N = No + ^l d\-g3/2 (^e-^^^--)) (14) 

In this intuitive approach, however, the Bose condensate was not described accurately. 
As it turns out, the problem can be solved exactly [5]: 

n{r)^-^\Mr)f + ^G{z,r) (15) 

1 — Zi A'' 

where ipo{r) is the ground-state wave function in the potential, and 

r.f ^ 1 ^ J exp[-(r/ro)nanh(e^/2)] ^ \ 

= I [l-exp(-2.£)]3/^ -exp[-(r/ro) ]| (16) 
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with 



hu/kT 

^A;r/27rma;2 (17) 

The exphcit occurence of ■0o('") shows that Bose condensation occurs in the ground state ipo 
of the potential. The zero-momentum state is irrelevant here. The total number of particles 
is 

N = + ^ / d^rG{z, r) (18) 

1 — Zi J 

which shows that the condensation is a continuous process, though it may appear to be 
abrupt, when N is so large that the first term can be neglected except near zi — 1. 
The Thomas- Fermi approximation is good when e ^ 1. In that case we have 

G{z^,r)^g,/,{z^e-^''^^^) (19) 

and therefore 

n(r) « ^|V;o(r)|2 ^ ^^^^^ (^le-'^^W) (20) 

which is similar to the naive formula, except for a better representation of the condensate. 
(The replacement of z by Zi is inconsequential.) The lesson is that a purely intuitive approach 
is not satisfactory, and we need a systematic method. 



Zi = 

e = 



3 Uniform Dilute Interacting Bose Gas 

The underlying idea of the Thomas- Fermi approach is to treat a nonuniform condensate as 
locally uniform, with a slowly varying density. I will first review the properties of a uniform 
Bose gas in the dilute limit, with interparticle interactions taken into account through a 
scattering length a > 0. 

The annihilation operator of a particle of momentum hk satisfies the commutation 
relation 

[ofe, aj.,] = 5kk' (21) 
We make a Bogolubov transformation to quasiparticle operators rjk'- 

ak = XkTik - VkV-k (k 7^ 0) (22) 
and require that the transformation be canonical, i. e. 

[Vk,Vk'] = ^kk' (23) 

This leads to the condition 
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which can be satisfied by putting 



Xk = cosh cTfc 

yk = sinhdfc (25) 



This a convenient parametrization, because interesting quantities find simple expressions: 

Pk = {alttk) = ^[cosh(2(Tfe) - 1] 

Afc = -(afca_fc) = ^sinh(2(Tfc) (26) 
where pk measures the depletion of the unperturbed condensate: 

No^N-Y,Pk (27) 

and Afc, is a measure of off-diagonal long range order. 

In the Bogolubov method, the annihilation operator for the zero-momentum state ao is 
equated to the c-number y/N. Explicit solution of the problem gives 

where p, is the chemical potential: 



Ai^ah^n /_ 32 na^\ 
m \ 3 V TT / 

with n the particle density. Note that this cannot be continued to negative a; apparently, 
new physics arises when the scattering length turns negative. The excitation energy of a 
quasiparticle of momentum p is given by 



.2 \2 



4 Quasiparticle Field Operator 

In the uniform case, the field operator ^'(r) can be put in the form 

^'(r) = ao + V(r) (31) 

with 

^^(0=0-1/2 5^ afce*-- (32) 

where Q, is the spatial volume. We have 

[^(r),^t(r')] = 5(r-r') (33) 
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since oq is treated as a c-number. We can introduce a quasiparticle field operator: 

Note that the relation between ip and ^ is non-local: 

V'(r) = / A[X(x - y)e(y) - y*(x - y)^t(y)] (35) 

where 

y(x-y) = ^ |/fee*-(--^) (36) 

For a non-uniform Bose gas, we write 

*(r) = 0(r) +^(r) (37) 
where 0(r) is a c-number function, such that 

(^(r)) = (38) 

where () denotes ground state expectation value. We transform to quasiparticle operators 
by putting 

V^(r) = J d'y[X{^, y)e(y) - F*(x, y)^\y)] (39) 

The requirement 

[e(r),et(r')]=5(r-rO (40) 

leads to the condition 

I dh[X{^, z)X(z, y) - F(x, z)r (z, y)] = 5(x - y) (41) 

The fulfillment of this condition in a simple fashion will guide our formulation of the Thomas- 
Fermi approximation. 

5 Wigner Representation 

In quantum mechanics, the Wigner distribution associated with a wave funcition ip{r) is 
defined by 

Pw(R, P) = / d^ri>*{R + r/2)V^(R - r/2)e'^-^/'' (42) 

That is, we take the off-diagonal density at two different points in space, and Fourier analyze 
with respect to the relative separation. This is illustrated in Fig.3. 
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Figure 3: To get the Wigner distribution, Fourier analyze with respect to relative coordinate. 



The Wigner distribution is not positive-definite, and hence not a probability; but it acts as 
a quasi-distribution function in phase space: 

(V^,pV^) ^ j dhr{r)\v^{v) = j'^^^VpwiR.p) (43) 



For a function X(x, y) that depends on two spatial points, we define its Wigner transform 

as 



.(R, P) = / dhX{R + r/2, R - r/2)e^P-'' (44) 
with the inverse transform 



If C (x, y) has the form 

C7(x,y) = J dhA{^,z)B{z,y) (46) 
then its Wigner transform takes the form 



1 ^ / 

C^(R,p) = Aw(R,p)Bw(R,p) + -J2y 



dAw dBw dBw dAw 
dRj dpj dRj dpj 



iJ^/ d^Awd^Bw d^Bwd^Aw d^Aw d^Bw ^ 
8^A^^? ^P? ^ dp'i dRjdpj dRjdp] 



+ ■■■ (47) 

The second term is the classical Poisson bracket {Aw, Bw}pb- It and the subsequent terms 
all depend on spatial derivatives, and would be small if the system is nearly unform. Thus 
our version of Thomas- Fermi approximation consists of keeping only the first term. Errors 
incurred can be estimated by calculating the next non-vanishing term. 
In terms of the Wigner transform, we can write 

I dhX{x, z)X{z, y) = Xw{R, p)Xw(R, p) + Yi^Xw, Xw}pb + ■■■ (48) 
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where the second terms vanishes identically. Thus, for such an integral, errors incurred in 
using the Thomas- Fermi approximation starts with subsequent terms. The condition on X 
and Y therefore reads 

X^(R,p)-Fi^(R,p)«l (49) 

and is solved by setting 

Xi^(R, p) = cosh(7(R, p) 

Yw(R,p) = sinh(7(R,p) (50) 

This make the problem very similar to the uniform case. 

At zero temperature, the criterion for the validity of the Thomas-Fermi approximation is 

huj/ii < 1 (51) 

where huj is the characteristic energy of the external potenial, and fi is the chemical potential. 
For the dilute interacting Bose gas, is of order of the scattering length. Thus, the Thomas- 
Fermi approximation can be used only when there are interparticle interactions. 



6 Variational Calculation 

We study the system defined by the Hamiltonian H, with 

H -nN = J d^x^^/i* + i y d^xd^y^^^ (52) 
where V{x) is the interparticle potential, and 

h^-—V^ + V,^{x)-pL (53) 
2m 

with V^xt(x) the external potential. The ground state free energy is 

F^{H- i^N) (54) 

where () means expectation value with respect to the ground state oi H — fxN . As mentioned 

before, we displace the field by writing = (p + where is a c-number function, such 
that {ip) = 0. We assume a trial form for the ground state, so that {F) has the same form 
as in mean-field theory, i.e., we can put 

(V'^(y)V^(x)V(x)V(y)) = A*(y, x) A(x, y) + p(y, x)p(x, y) + p(y, y)p(x, x) (55) 

where 

p(x,y) = (^^(x)^(y)) 

A(x,y) = -(^(x)V;(y)) (56) 

The ground state free energy F[(f),p,A; /i] is a functional of 0, p, and A, and also depends 
on /i as a parameter. The requirement {ip) = means that there are no terms in F linear in 
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Although we do not need the trial state explicitly, it can be explicitly constructed if 
desired. One can show that the wave functional of this state is of Gaussian form Thus, 
we have a true variational problem. 

We rewrite the functions p and A in p, A;/i] in terms of their Wigner transforms, 
and implement our version of the Thomas-Fermi approximation, as explained before. We 
transform to quasiparticle field operators, and find that, as in the uniform case, pw and A 
are parametrized by a single function: 



Ph'(R, p) 
Aiy(R,p) 



-[cosh(2a(R,p))-l] 
^sinh(2a(R,p)) 



(57) 



The Free energy reduces to the form F = J d^Rf(R.). We obtain equations for a and by 
minimizing F. The equation for is a modified Gross-Pitaevskii or non-linear Schrodinger 
equation (NLSE): 



■— + V;.t(r) + f/(r) - /i + v{0)<j)\r) 
2m 







(5^ 



where U{r) is a self-consistent potential that depends on a. It is unimportant for low 
densities. 



7 Dilute Interacting Gas in Harmonic Trap 

I will just quote some results for a dilute gas in a harmonic trap. The external potential is 

v«=%;r;V (59) 



2 \L 

For particles of mass m, 



^h/muo (60) 



which is the extend of the ground state wave function. For the interparticle interation, we 
use a pseudopotential 

d{v)—r (61) 

m or 

The sole effect of the differential operator above is the removal of a divergence in the ground 
state energy. The three important lengths in the problem are 

L (Extend of ground state wave function) 
a (Scattering length) 
i?o (Extend of condensate) (62) 

They are illustrated in Fig.4. 
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condensate 

ground state 



Figure 4: Length scales in atomic trap. Groundstate wave function is il^o{r). Condensate 
wave function is (f)o{r). 



For low densities, the non-linear coupled equations for a and (f) are solved by iteration, 
and one iteration suffices. The chemical potential is found to be 



\2lJ 60 V l y 



The requirement Tiu)/ fi <^ 1 means 



aN 



< 1 



The extend of the condensate is given by 



T ~ \ L ) 



(63) 



(64) 



(65) 



For the method to be valid we must have Rq » L. By neglecting the term V^^ in the 
NLSE, we find 

pa r / ^ \ 21 



STraL^ 



l-(- 
Rq 



(66) 



In Fig. 5 we show the shape of the condensate and estimated errors. Fig. 5(a) shows 0(r) 
as a function of r in units of L, for N=10^, and 10^. 

Fig. 5(b) shows the errors arising from the neglect of V^^. This is "trivial," as it can be 
corrected through numerical computation. 

Fig. 5(c) shows the errors incurred due to the Thomas- Fermi approximation, and are 
intrinsic to the method. They are small except at the edge of the condensate. 



10 



10 particles 

10^ particles 



4 6 
Distance (in units of L) 
1=10^* cm a = 5x10"^ cm 




10 



Figure 5: (a) Condensate wave functions for = 10^ and 10^; (b) Error incurred in neglect- 
ing kinetic term in NLSE; (c) Error incurred in Thomas- Fermi approximation. Length scale 
on horizontal axis is in units of L, the extend of the groundstate wave function. Calculations 
are done for L = 10~^ cm, scattering length=5 x 10"'' cm. 
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8 Quasiparticle Excitation 



The local excitation energy should be measured from the chemical potential: 



€p{r) = fx + 



\ [2m 




(67) 



where 



fJ'esir) = /i - Ve^t{r) 



(68) 



It describes a phonon with a position-dependent sound velocity. The excitation energy 
density of the system is given by 



where ej is the energy of the ith excited state. In the spirit of the Thomas-Fermi approxi- 
mation, we take 



The results for = 10^ and 10^ are shown in Fig. 6 and Fig. 7, with comparison to the ideal 
gas. 

Further details can be found in 

This work was supported in part by funds provided by the U.S. Department of Energy 
under cooperative agreement # DE-FC02-94ER40818. 



(69) 




(70) 
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Density of States for 1 particles 

1=10"* cm a = 5x10^ cm 




5 1^ 10 15 20 

Energy (in units of co/2) 



Density of States for 1 particles 



1=10^^0111 a = 5x10^ cm 




Energy (in units of co/2) 
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